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Abstract 



This paper studies the possibiUty of using the survival function, to predict 
long term stability by extrapolation. The survival function is a function of 
the initial coordinates and is the number of turns a particle will survive for a 
given set of initial coordinates. To determine the difficulties in extrapolating 
the survival function, tracking studies were done to compute the survival 
function. The survival function was found to have two properties that may 
cause difficulties in extrapolating the survival function. One is the existence 
of rapid oscillations, and the second is the existence of plateaus. It was 
found that it appears possible to extrapolate the survival function to estimate 
long term stability by taking the two difficulties into account. A model 
is proposed which pictures the survival function to be a series of plateaus 
with rapid oscillations superimposed on the plateaus. The tracking studies 
give results for the widths of these plateaus and for the seperation between 
adjacent plateaus which can be used to extrapolate and estimate the location 
of plateaus that indicate survival for longer times than can be found by 
tracking. 



Chapter 1 
Introduction 



This paper studies the possibihty of using the survival function, to predict 
long term stability by extrapolation |jl], ||, |], ^, H, ^. The survival func- 
tion is a function of the initial coordinates and is the number of turns a 
particle will survive for a given set of initial coordinates. To determine the 
difficulties in extrapolating the survival function, tracking studies were done 
to compute the survival function. The survival function was found to have 
two properties that may cause difficulties in extrapolating the survival turns 
function. One is the existence of rapid oscillations, and the second is the 
existence of plateaus.lt was found that it appears possible to extrapolate the 
survival function to estimate long term stability by taking the two difficulties 
into account. 

A model is proposed which pictures the survival function to be a series of 
plateaus with rapid oscillations superimposed on the plateaus. The tracking 
studies give results for the widths of these plateaus and for the seperation 
between adjacent plateaus which can be used to extrapolate and estimate 
the location of plateaus that indicate survival for longer times than can be 
found by tracking. 
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Chapter 2 

The survival function, jtms- 



For a given set of initial coordinate, xo,pxo,yQ,Pyo, one can find by tracking 
the survival time in turns, which is the number of turns the particle will 
survive before becoming unstable and which will be denoted by jtms- This 
determines the function jtrns{xo,Pxo, yo,Pyo) which will be called the survival 
function 0, |ll|]. If one limits the tracking to 1.0 10^ turns or less, one 



can find jtms for those xo,pxo,yo,Pyo for which jtms is less than or equal to 
1.0 10^ 

A tracking study was done of particle motion with no rf present, using 
an older version of the RHIC lattice. Random and systematic field error 
multipoles are present up to order 10. The particle momemtum is dp/p = 0. 
As the first case studied, the initial coordinates xo,pxo,yQ,PyQ are chosen 
along the direction in phase space given by pxo = 0,Pyo = and exo = e^o, 
where e^o, ^yo are the linear emittances in the absence of the error multipoles. 
Along this direction in phase space , jtms may be considered to be a function 
of Xq. For a given initial coordinate, Xq, one can find by tracking the survival 
time in turns, which is the number of turns the particle will survive before 
becoming unstable and will be denoted by jtms- This determines the survival 



function ||, [1^, |TI| jtms{xo)- If one limits the tracking to 1 10 turns or 
less, one can find jtms for those xq for which jtms is less than or equal to 
1 10^. For dp/p = 0, it is assumed that on the closed orbit x = and xq is 
also the initial betatron oscillation amplitude. 

The tracking was first done using xq which are seperated by dxo = .1 
mm. The results for jtms versus Xq are shown in Fig. ^]l|. The results in this 
figure may be looked at as the results of a search in xq starting at large Xq 
and decreasing xq in steps of dxo = .1mm . The figure shows an apparent 
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jtrns vs xO , dxO=.1mm, epxO=epyO 
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Figure 2.1: jtms versus Xq. dp/p = 0, ^^.q = 0,Pyo = 0,(^xo = e^o direction, 
dxo=.lmm. In the figure, jtrns, xO , dxO, epxO, epyO represent jtms, Xq, dxo 




stability limit for le6 turns of Usi =15.2 mm. 



2.1 Rapid oscillations in the survival function 



Fig. 2]1] shows rapid oscillations in jtms with xq. Large changes in jtms 
occur when xq changes by .1 mm. The oscillations extend over about .3 
mm. Tracking studies show that the oscillations become more rapid when 
the search interval, dxo, is decreased. This is indicated in Fig. p]2| where dxo 



is decreased to .05 mm and in Fig. |2.3| where dxQ is decreased to .025 mm. 
Results for different search intervals, dxo, are shown in Table |2.1| . dxo is 
decreased from .1mm to .0001 mm . The wavelength of the oscilations , Axq, 
as measured from the stability limit for le6 turns, Ugi, to the location of the 
first peak in jtms decreases from about .3 mm to .0002 mm. Also listed for 
each value of dxQ are the apparent stability limit for le6 turns, Usi, and jtms 
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dxQ (mm) 


Axq (mm) 


Usi (mm) 


jtrns Usl + dXQ 


.1 


.3 


15.2 


100000 


.01 


.06 


15.27 


150000 


.001 


.003 


15,270 


82605 


.0001 


.0002 


15.2793 


60693 



Table 2.1: Results for different search intervals, dxo- dxo is decreased from 
.1mm to .0001 mm . Axq is the wavelength of the oscilations , as measured 
from the stability limit for le6 turns, Ugi, to the location of the first peak in 
jtrns ■ Also listed for each value of dxQ are the apparent stability limit for 
le6 turns, Usi, and jtms at Ugi + dxQ. 



at Usl + dxQ. 

Table ^]l| shows that the wavelength of the oscilations , Aa;o is roughly 
proportional to the size of the search interval, dxQ. The value of jtms at 
Usl + dxo shows that near Usi, jtms changes appreciably in the small change 
in xq given by dxQ. Tracking results show that this seems to hold even at 
extremely small dx^.The computed results appear to indicate that jtmsi^o) 
is not a continuous function of xq. For a continuous function of xq, one can 
find a small enough interval in xq such that the differnce between the values 
of the function for any two Xq in that interval is very small. This does not 
appear to be true for jtrns{xo)- 

The existence of the rapid oscillations in the survival function, jtmsixo), 
would seem to make it difficult to extrapolate to find those xq that sur- 
vive for more than le6 turns. However one could view the survival function 
shown in Fig. ^]l| as being made up of the rapid oscillations superimposed 
on a smoother, more slowly varying function which could be used for the 
extrapolation. This is discussed further in section 2.3. 



2.2 Plateaus in the survival function 



Looking at Fig. |2.1| , one can make out plateaus in the survival function, 
jtmsi^o)- The plateaus are regions where jtms oscillates rapidly around an 
almost constant value of jtms- The plateaus can be seen somewhat more 
clearly if one reduces the search interval dxo, as shown in Fig. where dxo 



is decreased to dxo = .05mm. One can make out 4 plateaus located at about 
jtrns=^-5e5, 2e4, 3500, 400 turns. Possible plateaus with jtms less than 100 
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jtrns vs xO , dx0=.05mm, epxO=epyO 
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Figure 2.2: jtms versus xq. dp/p = 0, ^^.q = 0,Pyo = 0,(^xo = e^o direction, 
(ia;o=.05mm. In the figure, jtrns, xO , dxO, epxO, epyO represent jtms, Xq, dxo 

turns are being ignored. It will be seen that the width of the plateaus do not 
depend on the search interval, dxQ. This is also true of the location of the 
plateaus in jtms- This is shown in Fig. |2.^ where dxQ=.025 mm. 

The existence of the plateaus in the survival function, jtmsixo), would 
seem to make it difficult to extrapolate to find those Xq that survive for more 
than le6 turns. If one does the extrapolation using points which are close to 
the stability limit for le6 turns, Ugi, one may get incorrect results as these 
points may lie on one of the plateaus. 

2.3 Extrapolation of the survival function 

The data given above leads to a model of the survival functioon, which 
pictures it as sequence of plateaus. Within the plateaus, jtms oscillates about 
some constant vaue of jtms which will be called the level of the plateau. The 
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jtrns vs xO , dxO=. 025mm, epxO=epyO 
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Figure 2.3: jtms versus xq. dp/p = 0, ^^.q = 0,Pyo = 0,(^xo = e^o direction, 
da;o=. 025mm. In the figure, jtrns, xO , dxO, epxO, epyO represent jtms, Xq, 

dXo ,6x0, CyO- 

existence of the plateaus makes extrapolation of the survival function appear 
difficult. The last plateau that was measured has a level of about 1.5e5 turns 
and a width of about 1.2 mm. An interesting question is what is the level of 
the next plateau, at lower xq, after the last plateau that was measured. In 
the following this plateau will be referred to as the 'next plateau'. 

Some help with the extrapolation is also provided by plotting 1/ Log{jtrns) 
against xq as shown in Fig. p.4| , where the search interval is (ixo=.lmm and 
points with jtrns less than 400 turns are omitted. It will be seen below 
that the seperation between adjacent plateaus does not vary greatly when 
measured as the change in 1/ Log (jtms)- 

To help locate the 'next plateau', long runs were done starting with 
xo=15.2mm, and decreasing xq In steps of .1mm. In order to detect the 
beginning of the 'next plateau', the runs have to be long enough not to be 
confused by the rapid oscillations in jtrns that occur within each plateau. 
Runs of length 2e7 turns were chosen, and these runs take about 10 days for 
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1/Log(jtrns) vs xO , dxO=.1 mm, epxO=epyO 
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Figure 2.4: l/Log{jtrns) versus xq. dp/p = 0, p^ = 0,Pyo = 0,6^ = e^o 
direction, (ixo=.lmm. In the figure, Log(jtrns), xO , dxO, epxO, epyO represent 

trns 



the RHIC lattice used. 

The resuhs are shown in Fig. |2.5| . One sees from this figure that the 'next 
plateau' appears to start at a;o=15.3 mm and ends at about 14.0 mm, and 
has a level of about jtrns=15e6 turns and a width of 1.2mm . The previous 
plateau goes from 16.5mm to 15.2 mm and has a width of 1.3mm and a level 
of jtrns=1.5e5 turns. The width of the plateau is measured here from the 
end of one plateau in xO to the end of the adjacent plateau and includes the 
transition region where the points move from one plateau to the next. The 
width of the 'next plateau' is difficult to measure, as the adjacent plateau at 
still lower xq is estimated to have a level of 3e9 turns, and cannot be found 
by tracking. One can see that the width is larger than 1.2mm. The xq at 14.0 
mm survived more than 89.9e6 turns. The data given above for the 'next 
plateau' is somewhat in error but it seems better to use it than to throw 
away this information. 

The results are also shown as a l/Log{jtrns) versus Xq plot in Fig. |2^ . 
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jtrns vs xO , dxO=.1mm, epxO=epyO 
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Figure 2.5: jtms versus xq including points with jtms up to 4e7. dp/p = 0, 
Pxo = 0,Pyo = 0,ea;o = eyo direction, dxo=.lmm. In the figure, jtrns, xO , 
dxO, epxO, epyO represent jtms, xo, dxo ,6^, eyo- 

One can see three plateaus near the stability boundary at the levels of about 
jtrns= 15e6, 1.5e5,and 2e4 turns. The properties of the plateaus are summa- 
rized in Table |2.2| . The seperation between the plateaus in the 1/ Log{jtrns) 
plot are given by .054 and .039, so that the seperation between the plateaus 
in the l/Log{jtrns) are not too different. The widths of the plateaus that 
were measured were 1.2+, 1.3 and 1.7 mm. 

The plateau model will now be used to extrapolate and investigate long 
term stability in RHIC. In RHIC, jfrns=le9 turns coresponds to a survival 



plateau level, jtms 


15e6 


1.5e5 


2e4 


plateau level,!/ Log (jtms) 


.139 


.193 


.232 


plateal seperation in 1/ Log{jtms) 


.054 


.039 




plateau width, Axq (mm) 


1.2+ 


1.3 


1.7 



Table 2.2: Plateau parameters for the epxO=epyO direction. 
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1/Log(jtrns) vs xO , dxO=.1mm, epxO=epyO 
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Figure 2.6: 1 / Log{jtrns) versus Xq showing the next plateau and the last mea- 
sureed plateau, dp/p = 0, p^o = 0,Pyo = 0,6^0 = e^o direction, dxo~.lmm. 
In the figure, l/Log( jtrns), xO , dxO, epxO, epyO represent 1/ Log{jtrns), xq, 

dXo ,6x0, ^yO- 

time of about 3.5 hours. The data given above will be used to extrapolate 
to find the plateau whose level is greater than or equal to jtrns=^^^ turns. 
This plateau will be called the leQ plateau. The existence of the plateaus 
indicates that there are limits on the accuracy one can hope to achieve by 
extrapolation. The two parameters one needs to extrapolate the survival 
function are the plateau width and the plateau level seperation. One cannot 
be certain what these parameters will be at jtrns of the order of le9 turns. 
However, one can use the data found for these two parameters at jtms of the 
order of le6 turns , to make the best estimate for these parameters at jtms 
of the order of le9 turns. 

The plateau width and the plateau level seperation were studied for three 
different cases coresponding to three different directions in phase space. The 
results are summarized in chapter 5. The widths of the plateaus were found 
to be roughly constant when measured in terms of Xo=\/3xo{6xo + Cyo)]'^ with 
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Figure 2.7: 1/ Log{jtrns) versus Xq showing the plateaus found by extrap- 
olation, including the le9 plateau, dp/p = 0, p^o = 0,Pyo = 0,6^0 = ^yo 
direction, (ixo=.lnim. In the figure, l/Log( jtrns), xO , dxO, epxO, epyO 
represent 1/ Log{jtrns), Xq, dxo ,e^o, e^o- 



an average value AXo=2.00 mm. The plateau level seperation measured in 
l/Log{jtrns) varied from .054 to .019 with an average value of .033. It is sug- 
gested that in the extrapolation, the average value of these two parameters 
be used. For the case being considered here, it is assumed in the extrap- 
olation that the plateau widths in xq will be 1.4 mm corresponding to the 
average value of AXo=2.00 mm, and the plateau level seperation is .033 in 
l/Log{jtrns)- Note that in this case Xo=1.414 xq 

The results found using these assumptions are shown in Fig. p.7| , where 
two new plateaus are shown that were found by extrapolation. The plateau 
just below the next plateau in xO has the level of jfrns=3.3e9 turns and is 
the le9 plateau in this example. The le9 plateau goes from xo=13.9mm to 
Xo=12.6mm. The aperture for le9 turns may be taken as about 13.9mm, 
which may be compared to the 15.2 mm found using runs of le6 turns. This 
indicates a loss of about 1.3mm or about 9% due to the required survival 
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time of le9 turns. One may note that the level of the le9 plateau being 3.3e9 

turns, one may expect that some of the xq on this plateau will not survive 
leQ turns due to the oscillations in jtrns that will occur on this plateau. To 
be safer one could assume the aperture for leQ turns to be 12.6 mm, which 
is the beginning of the adjacent plateau with a level of 1.2el4 turns, giving 
a loss of 17%. 

The procedure used indicates that in this case the result is not sensi- 
tive to changes in the two assumptions made regarding the width and level 
seperation of the plateaus. If the plateau model continues to hold for the 
extrapolated plateaus, and the width and level seperation remain roughly 
1.3 mm and .034 respectively, then the result for the aperture for leQ turns 
will be about the same. 

The plateau model developed above avoids certain problems that arise in 
tracking studies. In trying to find the aperture for the survival time of le9 
turns, the more usual approach is to try find the Xq such that all smaller xq 
will survive for le9 turns. If while doing a tracking search starting from large 
Xo, one finds a Xq that survives for leQ turn, then one has to ask whether all 
smaller Xq will survive for le9 turns. This is a difficult question to answer. 
It is even possible, that there are are always smaller xq that will not survive 
le9 turns, although these xq may become very scarce at smaller xq. In the 
plateau model, in trying to find the aperture for the survival time of le9 
turns, the approach is to try to find the plateau whose level is greater than 
le9 turns. This is a better defined target, and the expectation is that this 
plateau will indicate a region of xq where most xq will survive le9 turns. 
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Survival function along the 
eyQ = direction 
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Figure 3.1: jtms versus xo ■ dp/p = 0, pm = 0,Pyo = 0, e^o = direction, 
dxo=.lnim. In the figure, jtrns, xO , dxO, epxO, epyO represent jtms, xq, dxo 
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jtrns vs xO , dx0=.05 mm, epyO=0 
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Figure 3.2: jtms versus xq . dp/p = 0, p^o = 0,Pyo = 0, e^o = direction, 
(ia;o=.05mm. In the figure, jtrns, xO , dxO, epxO, epyO represent jtms, Xq, dxo 

The tracking results presented above were all done along the exo = Cyo 
direction. Results will now be given for the yo = 0,Pyo = or e^o = 
direction. Tracking studies of other cases will test the consistancy of the 
plateau model. The direction in the space of Xo,Pxo,yo,Pyo is further defined 
by Pxo = 0. The particle motion is 4 dimensional because of the presence 
of skew multipoles. jtms may then be considered to be a function of Xq. 
Fig. shows jtrns as a function of Xq as found with tracking runs of le6 
turns. The apparent stability limit using le6 turns and dxQ=.lmm is 22.0 
mm. To make the plateaus more visible , one can reduce the search interval 
dxo to dxo=.05 mm. These results are shown in Fig. |3^ . One sees that there 
are two plateaus in the region shown with jtms greater than or equal to le4, 
whose levels are located at jtms= -SeG , and le4 turns .The oscillations on 
the plateau near the stability boundary, appear to be smaller than those seen 
in the e^o = Cj/o case. 

Runs of about 2e7 turns were done to find the shape of the 'next plateau' 
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jtrns vs xO , dxO=.1 mm, epyO=0 
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Figure 3.3: jtms versus xq including points with jtms up to 4e7. dp/p = 0, 
Pxo = OjPjyo = 0, e^o = direction, dxo=.l mm with dxo=.05 mm at larger 
xo- In the figure, jtrns, xO , dxO, epxO, epyO represent jtms, xq, dxo ,exo, ^yo- 



in the region where Xq is less than or equal to 21.0 mm. Enough tracking runs 
were done to determine the beginning , and the level of the 'next plateau'. 
The results are shown in Fig. ^ . Including these longer runs , one sees that 
the 'next plateau' begins at x0=21.0mm and the level of the 'next plateau' 
is about jt.rns=8e6 turns. Here, the begining of the plateau is the edge at 
larger xq and the end is the edge at lower xq. Th end of the 'next plateau' 
is some what difficult to determine, as the adjacent plateau at lower xq has 
a high level of about le9 turns. The end of the ' next plateau' was taken to 
be at a;o=19.0 mm. The adjacent plateau at higher xO is at the level of 1.5e6 
turns and with the width of 1.8 mm. These two plateaus are seperated by 
.017 in l/Log{jtrns) which is smaller than the .054 found in the exo=€yo case. 
The results are also shown as 1/ Log (jtms) versus xq plot in Fig. |3^ . 

The data found for these two plateaus will be used to extrapolate and 
find the le9 plateau in this case. In the extrapolation, the results found in 
chapter 5 for the average plateau width and level seperation will be used. In 
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1/Log(jtrns) vs xO , dxO=.1 mm, epyO=0 
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Figure 3.4: l/Log{jtrns) versus xq . dp/p = 0, p^ = 0,Pyo = 0, e^o = 
direction, dxo=.l mm. In the figure, Log( jtrns), xO , dxO, epxO, epyO 
represent 1/ Log{jtrns), xq, dxo ,eM, e^o- 



this case, this leads to the assumptions that each plateau is 2.0mm wide and 
a plateau level seperation of .033 in 1/ Log{jtrns) will be used here. 

The results are shown in Fig. |3.5| , where two extrapolated plateaus are 
shown with the plateau levels 8.5e8, and 4.5el2 . The aperture for le9 turns 
was taken to be 17.1mm, the beginning of the plateau with a level of 4.5el2 
turns. 17.1 mm is to be compared with the aperture of 21.9 mm found with 
runs of le6 turns. A loss of 4.8 mm or 22%. 
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1/Log(jtrns) vs xO , dxO=.1 mm, epyO=0 
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Figure 3.5: 1/ Log^jtms) versus Xq showing the plateaus found by extrap- 
olation, including the le9 plateau, dp/p = 0, p^o = Q,Pyo = ^,^xo = ^yo 
direction, (ia;o=.lnim. In the figure, l/Log( jtrns), xO , dxO, epxO, epyO 
represent 1/ Log{jtrns), Xq, dxo ,6^, e^o- 



Chapter 4 

Survival function along e^o = 
direction 
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Figure 4.1: jtms versus Xo .Xq = {Pxo/ PyoY^Vo- dp/p = 0, a;o=0, p^o = 
^:PyO = Oj ^1:0 = direction, dxo=.l mm and runs go up to 2e7 turns. In the 
figure, jtrns, XO , dXO, epxO, epyO represent jtms, ^o, dXo ,e^o, ^yo- 
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1/Log(jtrns) vs XO , dX0=.1mm, epxO=0 
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Figure 4.2: l/Log{jtrns) versus Xq . dp/p = 0, p^o = 0,Pyo = 0,6^ = 
direction, dXQ=.l mm. Xq = {bexO/beyO)-^yo. In the figure, 1/ Log( jtrns), 
XO , dXO, epxO, epyO represent 1/ Log{jtrns), Xq, dXo ,€^0, Cyo- 



Results will now be given for the xq = 0,pxo = or e^^o = direction. In 
order to be able to compare results with those of the above two cases, jtms 
will plotted against Xq = {Pxo/ Pyo)'^yo- Fig. ^]T] shows jtms as a function 
of Xq as found with tracking runs of about 2e7 turns. The apparent stability 
limit using le6 turns and (ixo=.lmm is Usi=24.9 mm. In Fig. 



one can 



make out two plateaus with levels larger than le4 turns. The levels of these 
two plateaus are at 1.7e6 turns and 1.5e7 turns. The end of the 1.5e7 plateau 
was taken as Xo=21.4 mm. The results are also shown as a 1/ Log{jtrns) plot 



in Fig. g^. 

The data found for these two plateaus will be used to extrapolate and 
find the le9 plateau in this case. In the extrapolation, the results found in 
chapter 5 for the average plateau width and level seperation will be used. In 
this case, this leads to the assumptions that each plateau is 2.0 mm wide and 
a plateau level seperation of .033 in 1/ Log{jtrns) will be used here. 

The results are shown in Fig. E^, where two extrapolated plateaus are 
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1/Log(jtrns) vs XO , dX0=.1mm, epxO=0 

0.4 I , , , , , , , , — 
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I ' ' ' ' ' ' ' ' ' 
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Figure 4.3: 1/ Log{jtrns) versus Xq . dp/p = 0, p^o = 0,Pyo = 0,6x0 = 
direction, dXQ=.l mm, Xq = {bexO/beyO)-^yo, showing extrapolated plateaus. 
In the figure, 1/ Log( jtrns), XO , dXO, epxO, epyO represent 1/ Log{jtrns), 

Xq, dXQ ,6x0, fj/O- 

shown with the plateau levels 3.3e9, and 7.7el3 turns . Most of the particles 
on the 3.3e9 plateau will be assumed to survive lc9 turns, and the aperture 
for le9 turns is then 21.4mm. 21.4mm is to be compared with the aperture 
of 24.9 mm found with runs of le6 turns. A loss of 3.5 mm or 17%. 



Chapter 5 

Extrapolation parameters for 
the plateau model 



direction in phase space 


^yO—^xO 


ej/o=0 


exQ=0 


plateau level, jtms 


15e6,1.5e5,2e4 


8e6,1.5e6,le4 


1.8e7,1.7e6,4e3 


plateau level, l/Log{jtrns) 


.139, .193, .232 


.145, .162 


.138,.160,.277 


platea level seperation in 1/ Log (jtms) 


.054, .039 


.017 


.022 


average plateau seperation 
in l/Log{jtrns) 


.033 






plateau width, AXq (mm) 


1.70, 1.83, 2.38 


2.0, 1.8 


2.2+, 2.0 


plateau width, AXq/Xq 


.075, .051, .106 


.095, .086 


.092, .083 


average plateau width, AXq 


2.00 mm 







Table 5.1: Plateau parameters for three directions in phase space. Xq = 



The extrapolation depends on two parameters, the width of the plateaus 
and the seperation between consecutive plateau levels as measured in terms 
of l/Log{jtrns)- The behaviour of these two parameters was studied by do- 
ing tracking runs for 3 different cases, corresponding to the three different 
directions in phase space, eyo=0, exo=eyo and £^^0=0. The two parameters 
were measured for these 3 cases using runs of about 2e7 turns. The results 
are summarized in Table |5.1| . Altogether, 7 plateaus were found and the two 
parameters for these 7 plateaus were measured. In Table |5ri| one sees that 
the plateau width, as measured as AXq, Xq = {Pxoi^xo + ^yo))'^, is relatively 
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constant with an average value of AXo=2.00 mm. The plateau level sepera- 
tion when measured in 1/ Log{jtrns) varies considerably with an average value 
of .033. 

Based on the above results, it is proposed that in extrapolating the sur- 
vival function, it is assumed that the extrapolated plateau widths are given 
by AXo=2.00 mm, and the plateau level seperations in 1/ Log{jtrns) a-re .033. 



Chapter 6 
Conclusions 



Tracking studies lead to a model of the survival function, which pictures it 
as sequence of plateaus. Within the plateaus, the survival time in turns, 
jtrns, oscillates about some constant vaue of jtms which will be called the 
level of the plateau. Studying the survival function along different directions 
in phase space, using an older version of the RHIC lattice , one finds that 
the width of the plateaus, A.Xq, Xq = {Pxoi^xo + ^yo))'^, remains roughly 
constant at about 2.00 mm. The seperation between the levels of adjacent 
plateaus has the same order of magnitude when measured in terms of the 
change in 1/ Log{jtrns), and has an average vaue of .033. Using, these results 
for the width of the plateaus and the seperation between plateau levels, one 
can extrapolate to estimate the location of the plateaus that correspond to 
longer survival times than can be found by tracking. For the case treated, it 
was found that a required survival time of le9 turns reduced the aperture by 
about 15% as compared to the aperture found by tracking using le6 turns. 

The plateau model also leads to new criteria to be used in tracking studies 
to find the aperture for particles to survive a given number of turns. In the 
plateau model, one finds the first plateau whose level is higher than the given 
number of turns, in order to find the aperture for the given number of turns. 
This is to be compared with often used method, where one does a search 
starting at large amplitudes until one finds an amplitude that survives the 
given number of turns. In the latter method one cannot be sure that a finer 
search would not find unstable runs at smaller amplitudes or how frequently 
these unstable runs will occur. In the plateau model, while there may be 
unstable runs at smaller amplitudes, there is the assumption that they will 
not occur frequently. 
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